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Abstract 

Supercurrent flow is studied in a structure that in the Ginzburg-Landau 
regime can be described in terms of an effective double barrier potential. In 
the limit of strongly reflecting barriers, the passage of Cooper pairs through 
such a structure may be viewed as a realization of resonant tunneling with 
a rigid wave function. For interbarrier distances smaller than do = t^HT) 
no current-carrying solutions exist. For distances between do and 2(io, four 
solutions exist. The two symmetric solutions obey a current-phase relation of 
sin(A(/?/2), while the two asymmetric solutions satisfy Ac/? = vr for all allowed 
values of the current. As the distance exceeds ndQ, a new group of four 
solutions appears, each contaning (n — 1) soliton-type oscillations between 
the barriers. We prove the inexistence of a continuous crossover between 
the physical solutions of the nonlinear Ginzburg-Landau equation and those 
of the corresponding linearized Schrodinger equation. We also show that 
under certain conditions a repulsive delta function barrier may quantitatively 
describe a SNS structure. We are thus able to predict that the critical current 
of a SNSNS structure vanishes as \/T^ — T, where is lower than the bulk 
critical temperature. 
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I. INTRODUCTION 



The superconducting state is characterized by the existence of long-range phase coherence 
in the electron system, and its characteristic nondissipative currents are associated to spatial 
distortions of the macroscopic phase. Superconductivity is however not the only instance of 
electron transport being phase- coherent over distances much larger than atomic length scales. 
Another important example is provided by mesoscopic normal transport, which began to be 
investigated about fifteen years ago and which quickly matured into an active branch of solid 
state physics.S Mesoscopic transport is realized at low enough temperatures, when electrons 
may preserve its phase coherence over long distances and thus undergo multiple coherent 
scattering by impurities or boundaries. Interest in this area stimulated experimental and 
theoretical research on a rich variety of mesoscopic phenomena in solids,! each of which is 
associated to a specific quantum interference process. Let us consider for a moment one of the 
most characteristic examples: the Aharonov-Bohm (A-B) effect. This is the effect by which 
the electronic properties of a thin ring depend periodically on the threaded magnetic flux. 
As a mesoscopic phenomenon, it was first observed in narrow normal cylinders! and later in 
small rings.! However, the A-B effect in solids had already been observed in superconducting 
rings not long after the discovery of the Josephson effect .H This observation was possible 
without the availability of modern nanotechnology because of the existence of long-range 
phase coherence in superconducting rings. Phase rigidity of the superconductor collective 
wave function is enforced by the spontaneous breaking of gauge symmetry. In contrast, the 
spatial coherence of the electron field in a normal system cannot rely on the existence of a 
phase transition and requires low temperatures and short length scales to make itself noticed. 
One may adopt the simple picture that the A-B effect in a superconducting ring relies on the 
phase coherence of the Cooper pair wave function, while its normal ring counterpart requires 
coherence of the single electron propagator. This consideration leads us to the question of 
whether, for any given electron interference phenomenon observed in a normal mesoscopic 
system, there may be of a Cooper pair analogue that could eventually be observed in a 
macroscopic superconductor. We have seen that there is an A-B effect for single electrons 
as well as for Cooper pairs. In the first case, the conductance (more generally, the current- 
voltage characteristics) is a function of the flux with a strong periodic component.! In the 
second case, it is the critical current (more generally, the current-phase relation) that depends 
periodically on the flux. The question is whether other quantum interference phenomena 
may also display this dual nature whereby they can be realized in the propagation of single 
electrons or Cooper pairs. A proper comparison of the two types of dynamics requires the 
resolution of the equations satisfied by the corresponding effective wave functions. In the 
case of normal electrons, and within the noninteracting approximation, one has to look at 
the wave function of Fermi electrons. In the superconducting case, the attention must be 
turned to the order parameter ip{r) oc< ipiipi >, if one wishes to obtain information on the 
global condensate behavior. 

Within the Ginzburg-Landau (GL) approximation, ^(r) must be an extremum of the 
free-energy functional 

F = J rfr{|(V - zA)^r - (1 - V{v)M' + A|^|V2}, (1) 

where V{r) is an effective potential and reduced units have been used.!!! The parameter A 
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is introduced here to eventually distinguish between the nonlinear GL case (A = 
linear case of a Schrodinger electron with reduced energy equal to one (A = 0). 
solutions of Eq. (|1|) must satisfy the equation 



1) and the 
Stationary 



(iV + Af^ - (1 - V{v))i) + AIV-I V = 0. 



(2) 



From the study of Eq. (H), we plan to investigate how Cooper pairs behave under an 
effective potential V{v) that, in the absence of nonlinearity, is known to yield a specific 
quantum interference phenomenon. 

Resonant tunneling (RT) is one of the best understood and perhaps simplest quantum 
interference phenomenona.i It can occur in one dimension without the requirement of a non- 
trivial topology. A quantum particle undergoes resonant tunneling when it has to traverse 
a double barrier structure. As a consequence of multiple inner reflection, the transmission 
probability is very sensitive to such parameters as the electron energy or the interbarrier 
distance, peaking near well-defined resonances. In this paper, we choose resonant tunneling 
as a case study in which to analyze the possibility of a doublefold phenomenology in quan- 
tum interference processes. For simplicity, we focus on the case of a quasi-one-dimensional 
superconductor and assume 5-function barriers, so that V{y) in Eq. (|^) is taken to be of the 



In the case of g large and A = 0, Eqs. (0) and (0) yield the essential phenomenology 
of RT for independent electrons. For this reason, we will pay especial attention to the case 
A = 1 with (7^1, although the crossover to moderate values of g in the nonlinear case will 
also be analysed. 

We must specify the transport properties we wish to compute. Ideally, one would like 
to compare identical device properties, such as the I-V curve, which largely characterizes a 
normal RT diode. However, we already encounter at this stage a major difference between 
the two types of systems we wish to study. Unlike in the case of normal transport, the 
I-V characteristics is not the relevant quantity to obtain information about the Cooper 
pair dynamics. The application of an external bias V creates a nonequilibrium between the 
populations of quasiparticles coming from each side of the structure. Only very indirectly 
does this change in the quasiparticle population affect the properties of the condensate, and 
the existing effect is very sensitive to the djTiamics of single quasiparticles,! which is very 
similar to that of normal electrons. In contrast, the population of Cooper pairs cannot be 
driven out of equilibrium, since they can only exist in a condensed state. One could at 
most create a difference between the chemical potentials of the condensates on each side 
of a structure. However, in the case of two barriers, one only expects a straightforward 
realization of the ac Josephson effect. 

The current carried by the condensate is driven by spatial variations of the macroscopic 
phase. For this reason, we focus on the study of the current-phase relation for configurations 
■ip{x) satisfying Eq. (Q) with V{x) given by (^. The current-phase relation is only meaningful 
for a superconducting structure; thus, we can hope at best to give a qualitative comparison 
between the two types of RT. One might still foster hopes that something resembling the 
scattering wave function of a single electron might still be obtained in the formal limit of 
A ^ 0. As we will see, however, this is not case. Because of the difference in boundary 
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conditions, the solutions for the superconducting order parameter obtained from Eq. (H) 
with A = 1 and the wave functions for Schrodinger scattering states with A = fall into 
qualitatively different classes of solutions that do not evolve continuosly into one another as 
A ^ 0. 

The physical explanation to the major differences between the normal and supercon- 
ducting versions of RT is to be found in the existence of strong correlations among Cooper 
pairs, which is caused by their bosonic nature and by their large mutual overlap. The 
picture of independent particles undergoing coherent multiple scattering, which often ap- 
plies for Schrodinger electrons in the normal state, completely breaks down in the case of 
Cooper pairs. The notion of many independent pairs which individually experience quan- 
tum interference and collectively combine to form a macroscopic wave function would at 
best be adequate for a condensate formed by noninteracting molecular pairs. However, such 
a scenario is not known to occur in nature. One may have superfluids made of strongly in- 
teracting bosons (case of '^He) or largely overlapping pairs ('^He) , and superconductors with 
large (BCS) or moderately short coherence length (case of cuprate superconductors). The 
strong correlation in conventional superconductors manifests itself through the nonlinear 
term in Eq. (H). As is characteristic of broken symmetry states, the resulting wave function 
is rigid and does not obey a superposition principle. We will see that the stiffness of the 
macroscopic wave function is the primary cause of the important differences between the 
Cooper pair and the individual electron scenarios. 

Let us now return to the A-B effect. We have already noted that it can be observed both 
in superconducting and normal (mesoscopic) rings. Being a direct consequence of electro- 
magnetic gauge invariance,0 the A-B effect is of a very fundamental nature. The transport 
properties of a thin ring must always display a periodic (although not necessarily always 
observable) dependence on the magnetic flux. In contrast, the sensitivity of the current 
to, e.g., the interbarrier distance (for a given electron energy and for a given normaliza- 
tion factor) in a RT structure does not stem from any fundamental symmetry. We will see 
that a rich dependence of the structure of the set of solutions on the interbarrier distance 
does certainly exist. We will also conclude however that, at large separations, an important 
physical quantity like the critical current becomes essentially independent of the interbarrier 
distance. 

The search for a Cooper pair analogue of resonant tunneling has led us to investigate the 
interesting properties of superconducting structures that can be described by a free-energy 
functional of the type ([^). As is proved in Appendix A, under certain conditions, a quasi-one- 
dimensional SNS structure without current concentration (S and N have the same width) can 
be quantitatively described by a (5-function in the effective GL equation. As a consequence, 
the model given by Eqs. (0) and (0) can yield a quantitative description of a SNSNS structure 
of uniform width. We have found that the set of solutions to the corresponding GL equation 
can display a very rich structure. The main experimental prediction is that, in a SNSNS 
structure, the critical current vanishes with a law Ic{T) ~ ^Jt^ — T, which differs markedly 
from the {T(. — TY behavior of the SNS case.0 In addition, we find T'^ < T^, i.e., for a SNSNS 
system, there is a depression (with respect to the SNS case) of the temperature at which 
the critical current vanishes. 

This paper is arranged as follows. After presenting the model in section II, an analytical 
study of the solutions to the nonlinear Eq. (^) is presented in section III. In section IV, we 
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prove the inexistence of a continuous crossover between the physical solutions of Eq. 
and those of its linear version. Some experimental predictions on temperature dependent 
transport properties are discussed in section V. Finally, in section VI, we present some 
concluding remarks and comment on how some of the qualitative conclusions of the present 
article can be extended to other interference phenomena. 



II. THE MODEL 

In this paper, we wish to analyse the current- carrying solutions of the nonlinear GL 
equation with the effective potential given by (^. We may factorize the order parameter 
ip = Re^"^ and write a = A — Vy^ for the superfluid velocity. The gauge-invariant electric 
current is then written j = —R^sl. Within the assumption of a sufficiently narrow wire 
(width much smaller than the coherence and penetration lengths) it is safe to neglect the 
dependence of j and ip on the transverse variables. In these conditions, our analysis reduces 
to the study of the solutions of the nonlinear differential equation 

^ + [1 - gS{x - a) - g5{x + a)]R - ^ - Ai?^ = 0. (4) 

for arbitrary values of g (hereafter, A = 1, unless otherwise stated). In a quasi-one- 
dimensional superconductor, we may choose A = and write for the phase 

where the current density j is a conserved number, the total current through a lead of 
cross-section A being I = jA. For the superconducting order parameter we are interested 
in solutions satisfying the boundary conditions 



R — Roo 

(p{x) = qx ± Alp/2 



for X ±oo, (6) 



Roo being the biggest solution (i.e., that with the lowest free energy) of the equation 

R^-R^+ f = 0. (7) 

In Eq. (P), a nonzero value of g = j/R"^ accounts for a linear variation of the phase at 
infinity,t3 and the possibility of a phase offset 

/OO POO ]^ ]^ 

yi,)-,yi.^,,Jj—^-„, (8) 

has been introduced. As in the single barrier case,0 A(f will turn out to be a most convenient 
parameter to classify the solutions of Eqs. (^ and (|). 

One may wonder whether, apart from describing the condensate analogue of resonant 
tunneling, the double delta barrier model may quantitatively describe a specific physical 
system. Fortunately, the answer is yes. It is proved in Appendix A that, within the GL 
approximation, a 5-function effective barrier can serve as a quantitative model for a SNS 
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structure. More specifically, we show that, if is the coherence length of the normal 
metal, the effect of a normal segment of length L inserted in a superconducting wire can be 
rigorously modeled by a 5-function of strength 
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SAT 



(9) 



provided that L <^ <^ ^(T), where ^(T) is the temperature dependent coherence length. 
Since this effective delta barrier yields the correct matching properties of a realistic SNS 
system,lli we may conclude with confidence that a model with a double (5-function barrier 
will correctly describe a SNSNS system.0 In addition, the fact that both the barrier strength 
g and the effective length a are functions of temperature for a given physical structure (note 
that a = d/2C,{T), where d is the physical distance between the normal islands) permits us 
to use temperature as a convenient driving parameter to tune g{T) and a(T). This property 
has important experimental consequences that will be discussed in section IV. 

Regarding the temperature dependence of the effective parameters, it is interesting to 
note that, as T — Tc, one has g{T) oo and a{T) — > 0. Therefore, the Josephson limit 
{g ^ 1) can always be explored by driving T sufficiently close to T^. In the Josephson 
regime, the currents are necessarily much smaller than the bulk critical current. As a 
consequence, the spatial variation of the phase can be safely neglected for many purposes 
and the phase offset can be identified with the the conventional "phase difference" between 
the two superconducting terminals, A^p ~ yji — Lp2- In Refs.0lil it was shown that, for large 
g, the critical current in the presence of a ^-barrier is 



l/2g. 



(10) 



Since the bulk {g 
relation 



0) critical current is, in these units, jb = 0.385, one can derive the 



g = 1.30-ijb/i'^) 



(11) 



valid for large g. Going beyond the 5-barrier model, one can show that the critical current 
for a SNS system of arbitrary normal length L is0 



— smh( — , 



(12) 



provided that the resulting jc is much smaller than j^. One may readily note that the result 
(|l^) is consistent with Eqs. (|9|-p!oD in the limit L <^ 

III. EXACT SOLUTIONS 
A. General properties 



In this section, we study the solutions of Eq. (^ satisfying the boundary conditions 
]). There is a mechanical analogue that helps to understand some general properties of its 
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solutions.0 Eq. may be viewed as the force equation for a classical particle of unit mass 
with position R at time x moving under a potential 

v{R) = -XR^A + R^/2 + f/2R^. (13) 

In the GL case (A = 1), configurations represented by points A and B of Fig. la satisfy 
the asymptotic condition dR/dx = 0, with B energetically more favourable. The picture is 
that of a classical particle which, at remote times {x — oo), stays in point B. At negative 
times, it begins to roll down and, after receiving two kicks at times x = ±a, it returns 
asymptotically (x —>■ oo) to point B. At each pulse, the change in the mechanical energy 
e = i?'^/2 + v{R) is determined by the matching condition 

R'{x+) - R'{x-) = gR{x), (14) 

at points x = ±a. Fig. lb shows one particular solution R{x) that exactly correlates with 
the mechanical analogue schematically depicted in Fig. la. 

The energy before and after the two kicks is £ = v{Roo)- Eq. (p!4D indicates that the 
effect of barriers is that of making the velocity R'{x) more positive. If the system is to return 
to point B at X ^ oo, it is not difficult to see that the conditions feo < (where 6eo is the 
change in mechanical energy at point x = —a) and R'{—a~) < must be satisfied, and, for 
analogous reasons, R' (a-) > 0. These two constraints force the solutions to be ot the fotnH 

R^{x) = Z + U tanh^(A;(x — x_|_(_))) x > a (x < —a) 

62 - ei 



R^{x) = ei + (e2 - ei) sn^ ( a\ ^^^—^{x + a) + Fc 



63 - ei 



x\ < a (15) 



where Z is the smallest root of Z{Z — 2)^ = 8j^, U = 1 — 3Z/2, k = yjU/2, and sn(x|m) is a 
Jacobi elliptic function.0 Defining the function y{x) = R^{x) and the parameters dq = y{—a) 
and u'q = y\—a~) (both of them functions of x_), the matching condition (p^) at x = —a 
determines ei, 62, 63 as the roots of the polinomial 

P{s) = {s~ Z){s-l^ Z/2f + A5eos, (16) 

with ei < 62 < 63, and where the change in e is 

5e, = g{y'^ + gy,)/2. (17) 

Finally, 



Fn = F\ arcsin 



lyQ - ei 



V V ^2 - ei 



62 - ei 



63 - ei 



(18) 



where F{ip\m) is the incomplete elliptic integral of the first kind,tj and a in (|T5|) is defined 



as 



a = sgn{y'^ + 2gyo). (19) 
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In the limit (t/q + '2gyo) 0, it can be shown that both values of a in Eq. ([T5| ) lead to the 
same result for B?{x). 

In analogy with uq and y'^, we define yi = y{a) and y'l = y\a^). The second kick will take 
us asymptoticaly to the value R^o if it causes an energy change Ssi such that 5^1 + Seq = 0. 
This leads us to write the relationship 

y'l + gyi + y'^ + gyo = ^ (20) 

as the global matching condition. Eq. ( pO]) determines implicitly the parameter x_ to be 
introduced in Eq. (|15|). Through the integration of Eqs. (|-|^), each possible value of x_ 
uniquely determines one solution of Eq. @ satisfying the boundary conditions (H). Thus, X- 
is a parameter that completely characterizes a given physical solution. Our goal is therefore 
to solve numerically for all possible values of x_ satisfying Eq. ( ]20|) for a given value of j. 

A quantity of interest is the critical current, which we plot in Fig. 2 as a function of 
the (reduced) semi-distance between barriers for several values of the barrier strength. The 
result is given in units of the critical current for a single barrier with the same strength, 
Jf.{a, g) = Ic{a, g) / l'^\g). For separation distances much larger than ^(T) the critical cur- 
rent becomes identical to that of a single barrier, regardless of the value of g. This is the 
limit in which the two barriers are effectively decoupled. The decoupling at long distances 
is caused by the nonlinear term in Eq. (Q). In contrast, the dynamics of an electron obeying 
the linear Schrodinger equation is always sensitive to the presence of both barriers. We will 
see however that, in the nonlinear case, there are always solutions with high free energy 
that are sensitive to the presence of both barriers (in the sense that they cannot be viewed 
as simple combinations of single barrier solutions), but these energetic solutions become 
increasingly irrelevant at large separations. 

We observe in Fig. 2 that, for moderate values of g, there is a slight depression of the 
critical current with respect to the single barrier case, if the reduced semidistance is a ^ 1 
(a similar effect was noticed in Ref.lil for the case of weak links). This effect becomes more 
pronounced as g gets larger. For g = 10, the law l/2g already applies approximately for 
the critical current of single barrier and then Jc — 1/2 for a = 0, as expected from a single 
barrier with doubled strength. Deep in the Josephson regime {g = 20 or larger), the interval 
0<a<7r/2is practically depleted of solutions and the critical current is essentially zero. 
As will be seen later, this has noticeable experimental consequences, since the value of the 
reduced distance is temperature dependent. In the following two subsections, we analyse 
the main properties of the solutions of of Eqs. (§-|^) over the whole range of possible values 
of g. 



B. Josephson limit (g oo) 

In the limit of large scattering strength J_Gf—> oo), the presence of a single barrier is 
known to yield the ideal Josephson behaviorllHUlJ 

j = il/2g)sm{A^), (21) 

where the phase offset Aip can be identified with the conventional "phase difference" between 
the two weakly linked superconductors. For the double barrier case, the parameter 6e = —Seq 
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serves to characterize the continuous set of solutions in the regime of large g. For a given 
value of the current j, it is shown in Appendix B that Se takes values between and 
(1 — J'^)/4:, where J = 2gj is the current in units of the critical current for a single barrier 
structure. We also show in Appendix B that the matching condition (|20| ) can be rewritten 

as 

a = nG{6e), n = l,2,... 



where K{m) = F{^\m) is the complete elliptic integral of the first kind.ts For each integer 
n that meets the requirement (^) there are four solutions R{x), all of them corresponding 
to the same value of Se. The quantity 2G{6e) is the spatial period of the solution between 
barriers. Thus, for each solution of Eq. (^2]) that we may find for a given value of n in a 
structure of interbarrier distance 2a, we can always construct a solution for the structure of 
distance 2{a + G) whose index is n + 1 and which is identical to the previous one except for 
the presence of an extra oscillation. 

The structure of the solutions can be clearly appreciated in Fig. 3. In the Josephson 
limit, we see that there are no solutions for 2a < it [J^ = 0). For 2a > tt the critical current 
becomes nonzero and a group of four solutions appears, two symmetric and two asymmetric 
under the inversion x —>■ —x, the two asymmetric ones being the mirror image of each other. 
Their current-phase relation is shown in Fig. 4a for a = 1.74. The symmetric solutions form 
two branches that combine to yield the curve 

/ = /e(a,(?)sin(Av9/2), (23) 

which is the naive expectation for two Josephson junctions in series. In contrast, the two 
asymmetric solutions obey the law 

AyD = 71, for all /. (24) 

This peculiar current-phase relation can be simply understood if the two junctions are as- 
sumed to be sufficiently far apart. In this case, the central peak in R{x) between the two 
barriers becomes a long plateau and the two barriers behave as independent. For a given 
current, there are two possible values of the phase, Aip^ and A(y9_, satisfying Eq. (pi]). In 
the presence of the two barriers, these two phase values may combine in four different ways. 
In two cases, the phase difference is the same in both barriers. Then, A(p equals 2A(p^ or 
2A(/3_. This results in a sin(A(y9/2) law for the symmetric solutions. In the other two cases, 
a different phase change takes place in each barrier. But the total phase difference A(p is 
given by the sum of the two values, which is exactly 

Aif+ + Aif^ = 71. (25) 

This explains Eq. (|24D and the double vertical branch shown in Fig. 4a. The above set 
of considerations makes the results ( PBD and ( pi[ ) quite plausible, and even expected, for 
the case of large separations. The not so obvious result is that the same conclusions hold 
when the separation distance is comparable to C,{T) and the two barriers can in no way be 
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viewed as decoupled. This result can be proved rigourously from analytical considerations 
(see Appendix C). The smallness of the distance between barriers makes itself noticed only 
through a depression of the critical current with respect to the single barrier case, but not 
in the qualitative form of the /(A<^) curves. 

As the distance grows, the maximum current of the group of four solutions saturates to 
the critical value for a single barrier. As indicated in Fig. 3, a new group of four solutions 
appears when a > vr, and their current-phase relation is shown in Fig. 4b. Comparison with 
Fig. 4a reveals two features of the /(Ayj) relation: (i) for a not much larger than vr, the 
maximum current value of the second group of solutions has not yet reached the saturation 
value, and (ii) all the values of Ay9 are shifted by vr with respect to the first set of four 
solutions. Again, these features are easy to understand if the two barriers are assumed to 
be far apart. The new group of solutions resembles the first set of four, except in that a 
soliton (i.e., a spontaneous local depression of the gaplll) contributing an extra phase of vr 
has nucleated between the barriers. As the two barriers separate, the added depression of 
the order paramater, which may initially be viewed as the result of quantum interference 
oscillations, evolves into a well-defined, isolated soliton. Remarkably, the additional depres- 
sion in \ip\ contributes exactly vr to the phase, even when the two barriers are relatively close 
and the depression does not act as an isolated kink but rather as an additional oscillation. 
It must also be noted that, since j ^ 0, we have IV'I 7^ at all points. 

As the separation distance is made increasingly larger, the same pattern repeats itself. 
Each time the distance 2a exceeds an integer multiple of vr, a new group of four solutions 
appears, with a structure similar to that of the preceding set of four solutions except for an 
extra phase of vr. As the interbarrier distance continues to increase, new sets of solutions 
emerge periodically, always in groups of four. For a given distance, the groups of solutions 
that are most sensitive to the double barrier feature (as indicated by its not yet saturated 
maximum value) are those with a higher free energy resulting from a higher number of 
modulations in the order parameter. These very energetic solutions can be expected to be 
irrelevant in practice, except perhaps to account for fine details in the dynamic behavior. 
Therefore, it may be stated with reasonable confidence that the two barriers become effec- 
tively decoupled for practical purposes (also in regard to the dynamic behavior) when the 
maximum current of the second group of solutions has reached the saturation value. This 
happens approximately for 2a > 3tt. 

Computation of the critical current. Eq. (|2^) (see also Fig. 3) indicates that 6e is a 
unique function of a and n, so that the product in the r.h.s. hand side of Eq. ( |B1U| ) must 



be constant within a given group of four solutions. The highest possible value of J (which 
is always to be found in the n = 1 group) is obtained by imposing u^l — u) to take its 
maximum value of 1/4. Thus we conclude that the critical current Jc = 2gjc is given by the 
relation 

1/2 



K{'-^4^\- (26) 



1 - J2 / V 1 + a/1 - J2 



Eq. (pq) yields the curve Jdo.) for the critical current in the Josephson limit, and the result 
has been plotted in Fig. 2. 



10 



C. Intermediate and small values of g 



In Fig. 5 one can analyze the crossover between the regimes of large and small values of 
g for the same barrier distances we considered in the discussion of the Josephson limit (see 
Fig. 4). Let us comment the case of a = 1.74. As we depart from the limit of very large g, 
an extra branch appears of solutions with tt < A<y9 < 2tx and with a low critical value (see 
the graph in Fig. 5a for g = 20). As g decreases further, the low current solutions with Aip 
near 27r disappear and the resulting branch begins to shrink until it merges (for g ~ 0.2) 
with the solutions that formed the double vertical branch in the Josephson limit. For g ^ 0, 
the two sets of solutions with A(f ^ form a triple branch corresponding to three single 
solitons nucleated at points xo = and 



decreasing g, the branch of solutions with the smallest phase offset evolves towards the set 
of uniform solutions (those for which R{x) = i?oo and A(f = 0) that are characteristic of a 
perfect superconductor.lll 

As can be seen in Fig. 5b, the crossover presents similar characteristics for the case 
where a = 3.5. Like in the Josephson limit, the main differences for large g (with respect to 
previous shorter distance case) lie in the extra value of tt of Aip shown by the second set of 
solutions and in the presence of a four-fold branch with a saturated critical current value. 
As g decreases, the evolution of the four solutions with higher phase offset runs similar to 
that of their counterparts in Fig. 5a. New interesting features appear however for low values 
of g. At g = 0.5 and 0.2 it is clear that the two groups of four solutions begin to merge into 
a simpler pattern of single and double solitonic solutions (with A(p tending to vr and 27t as 
J — s> 0, respectively). Since double soliton solutions cannot exist in the transparent case, the 
corresponding (four-fold) branch begins to detach from the rest of curves and to decrease its 
critical current until it eventually disappears. For g —>■ the remaining branch is identical 
to its counterpart in the shorter distance case: three solutions of single solitons located in 
the same points as above and one uniform branch with zero phase offset. Incidentally, it can 
be shown that for a < (2fc)~^ ln(2 + y^) only the soliton at xq = survives. This is not the 
case however in any of the two distances considered here. 

An interesting consequence of studying the crossover between large and small values of 
g is that, among the curves shown in Fig. 5, we can recognize similarities with current- 
phase relations computed for other structures that, not being exactly SNSNS, share some 
common features. For example, Martin-Rodero et al.0 have performed a self-consistent, 
zero-temperature calculation of the current-phase relation in mesoscopic weak links, mod- 
elled by a tight binding chain linked to two broad Bethe lattices that act as superconducting 
reservoirs. In spite of the obvious differences between the two physical models, the similarity 
between some of the resulting curves is striking. For instance, the current-phase relation 
shown for the two longest chains in Fig. 2c of Ref.0 resemble some of the branches we obtain 
for g = 5 and 20 in the a = 1.74 case. It is interesting to note that the curves of Ref.0 that 
look alike correspond precisely to the case of strong internal reflection at the constriction. 




(27) 
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One concludes: (i) the nonlinear term of the GL equation has an effect very similar to that 
of self-consistency in a zero-temperature calculation; and (ii) the essential physics rests in 
the nonlinear effects taking place within the finite superconducting segment and the effec- 
tive scattering at its two ends, the physical details of the semi-infinite S leads being less 
important. Conclusion (i) is what one expects from inspecting the microscopic derivation of 
the GL formalism, and it has alreadv been noted in connection with the crossover for large 
to small g in the single barrier case.t3 



IV. INEXISTENCE OF A CONTINUOUS CROSSOVER BETWEEN THE 
LINEAR AND NONLINEAR PROBLEMS 



It has already been commented in the Introduction that there are qualitative differences 
in the physics described by Eqs. (2-4) in the linear (Schrodinger) and nonlinear (GL) cases. 
However, one might still think that, from a mathematical point of view, there could be a 
continuity (as a function of A) between the order parameter configurations satisfying the 
GL equation (A = 1) and the scattering wave function for a Schrodinger electron in the 
presence of potential (^ with A = and unit energy. Again, it is useful to consider the 
mechanical analog represented in Fig. 1. The wave function of a retarded scattering state 
is of the type ip ~ e*'^^ + re~'*^^, for x < —a, if the electron is coming from the left, so that 
the amplitude R = {ipl displays oscillations in the left asympotic region. The mechanical 
picture to introduce in Fig. 1 would be that of a particle which at remote times x —>■ —oo 
is oscillating between two return points with energy e < v{Roo), receives two sudden pulses 
at times x = =ba, and ends asymptotically in one of the stationary points (this corresponds 
to the uniform amplitude solution ~ te*^^ to be found in the x > a region). Thus we 
see that a major difference between the two cases lies in the asympotic behavior of R{x), 
which is always uniform for the superconducting order parameter while it may oscillate for a 
Schrodinger electron. The only remaining possibility to find a common mathematical wave 
function must be looked for in that combination of parameters which yields a transmission 
unity for the electron linear wave, since in this case R would be uniform on both sides of 
the structure. The equivalent mechanical particle would start from a stationary point to 
which it would return at late times. The problem is that, for an electron satisfying the 
linear Schrodinger equation, that stationary point has to be A (see Fig. 1), and not B, as is 
the case for the GL order parameter. Although a rigorous demonstration of this statement 
is possible, it suffices to note that, as A 0, the potential maximum at B moves toward 
infinite values of v and R (see B' in Fig. 1 for A = 0.8), effectively disappearing from the 
picture when A = 0. Thus the only stationary point that is available in the linear case to 
describe purely transmitted waves is precisely point A, qualitatively different from B even 
in the limit A ^ 0. We conclude that, as the nonlinear term in Eq. (^ is made to vanish 
(A ^ 0), there does not exist a continuous crossover between the physical solutions of the 
nonlinear Ginzburg-Landau equation and those of the linear Schrodinger equation. 



12 



V. EXPERIMENTAL PREDICTIONS 



The proximity between barriers is most noticeable when the low energy solutions are 
sensitive to it. Of course, the most dramatic effects can be seen for 2a < vr, when no solutions 
exist at all in the large g limit. This has interesting experimental consequences. We must 
first remember that a is the reduced semidistance in units of ^{T), i.e., a(T) = d/2C,{T). 
This means that, for a given structure with a fixed physical distance d, the reduced distance 
2a can be made arbitrarily small by driving the temperature sufficiently close to T^, since 
then ^(T) = C,q{1 — T/Tc)~^^^ — > oo (for a clean superconductor, = 0.74,^0) while for 
a dirty superconductor = 0.85^/T^, I and being the mean free path and the zero- 
temperature coherence length, respectively^) . If the double barrier structure is formed by 
two normal segments, we know from the analysis in section III that g{T) scales towards 
the Josephson limit. Therefore, there is a temperature < Tc above which d < 7r^(T). 
Since, at the same time, g{T) is very large, this has the consequence that no current-carrying 
superconducting solutions are allowed in the temperature interval < T < Tc. The result 
is that, for a SNSNS structure, there is a depression (with respect to the SNS case) in the 
critical temperature above which the critical current becomes zero. 

In Fig. 6, the critical current is plotted as a function of temperature for both a SNS and 
a SNSNS structure. The main feature of the resulting Ic{T), namely, the law with which it 
vanishes, is amenable to an analytical treatment. In section III, we derived Eq. (p6|), which 
determines Jc{T) for a given value of a(T). Since we are interested in finding out how the 
critical current vanishes, we may expand the r.h.s. of Eq. ( [26| ) for small values of Jc{T) and 
obtain 




^ "^^ + ^MT)]. (28) 



Tc 2 V 16 
The temperature T'^ for which Jc{T^) = is T^ = Tc — AT^, with 

AT, = (^^^ Tc. (29) 

Eq. (p9D is actually an upper bound to the value of Tc — T^, and it tends to the exact value 
when the condition 5'(T^') ^ 1 is consistently satisfied. 

After some simple algebra, we get for the critical current (in real units): 

Ic{T) = ll^\T':){^fl\^) f^V^V^' T (30) 



^37r2^ \ Tc 



where 



^'^ g{Tym' ^{T) ~ L^Qtic'^ ni'i' T^ ^'^^^ 

is the critical current of a SNS structure at temperature T, A being the cross area of the 
junction and k = X{T)/^{T) is the Ginzburg-Landau parameter of the superconductor. The 
square root behavior of Ic{T) for a SNSNS structure contrasts markedly with the (T^ — T)^ 
law of the SNS case.0 
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Simple estimates suggest that the predicted depression in the critical temperature should 
be measurable. For example, for a SNSNS structure made with superconducting Al (T^ = 
1.19 K and = 1-6 /im) with d = lO^o (which falls within the validity of GL regime) one 
would obtain a decrease of AT^ ~ 0.05 K. On the other hand, it is interesting to note that, 
if the condition d ^ is satisfied, there is a range of temperatures sufficiently far below Tc 
for which d > ^(T) and, thus, /c(T) ^ I^^^T). 

The sensitivity of the effective g to temperature can be exploited in other interesting ways. 
For example, by varying the temperature, one may drive a given SNS or SNSNS structure 
from the large to the small g regime. Consider a SNS structure made in a narrow wire 
of finite length at whose extremes we apply an external voltage. N is a normal metal with 
critical temperature Tc < Tc. When T gets close to Tc, g(T) becomes large, and so does ^(T). 
Both facts contribute to yield an ideal Josephson behavior (a large value of ^ (T) makes the 
length of the wire effectivelv shorter, thus facilitating the adiabatic response characteristic 
of the ac Josephson effectJi3). As an external voltage is applied, the current oscillates very 
rapidly and one observes a zero time-average, Jav = 0. As the temperature is lowered below 
Tc, the structure becomes of the type SSS with an effective g ^ I (we may choose S to be not 
very different from S). The situation is then close to that of a uniform superconductor and 
the system cannot respond adiabatically to the externally applied biaJilSii. Some type of 
resistive behavior has to be displayed, with the result that Jav 7^ 0. The net effect is that, 
by lowering the temperature, it is possible to drive a given SNS structure from adiabatic to 
resistive response. 

VI. CONCLUSIONS 

We have studied superconducting flow in structures which, in the Ginzburg-Landau 
regime, can be described by a double barrier effective potential. We have found that the 
critical current is depressed with respect to the single barrier case. In the limit of strongly 
reflecting barriers, Jc can become effectively zero for interbarrier distances d < do = tt^(T). 
As d exceeds do, four solutions (two symmetric and two antisymmetric) appear, as might be 
expected for two Josephson links in series. As d grows even larger, new solutions appear, 
always in groups of four. The supression of Jc for short distances has practical consequences, 
since do(T) oo as T ^ Tc, and the condition d < do is always reached by any device if 
T is sufficiently close to T^.. Thanks to this analysis, we have found that the critical current 
of a symmetric SNSNS structure is lower than that of a SNS structure with an equivalent 
N segment. More important, the law with which /c(T) vanishes is qualitatively different in 
each case (see Fig. 6): /c(T) decays as ^Jt[--T in a SNSNS structure, as opposed to the 
(Tc — TY behavior of the SNS case, with the inequality T^' < Tc being always satisfied. 

We have also shown that, because of the different boundary conditions, the solutions of 
the GL equation cannot evolve continuously towards the scattering solutions of the linear 
Schrodinger equation as the nonlinear term is formally made to vanish. 

We close this article by returning to the question that initially motivated it, namely, 
the possibility of finding a macroscopic quantum analog of a specific interference process 
such as the resonant tunneling of Schrodinger electrons. It was already commented in the 
Introduction that, given the important differences between the relevant transport properties 
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of the physical systems involved, the analogy would at best be qualitative. We have learned 
that a basic fact such as the dependence of transport behavior on the interbarrier distance - 
all other parameters being equal- effectively disappears in a superconducting structure with 
d 3> T^^{T) (see Fig. 2). Since ^(T) is the length scale needed for the nonlinear term in 
Eq. (0) to make itself noticed, we may adopt the view that the nonlinear term in the GL 
equation acts as a dephasing term that damps any interference effect requiring modulation 
of the wave function amplitude R (as would be caused, i.e., by the interference of waves 
traveling in different directions). This picture allows us to extrapolate our conclusion on 
the inequivalence of macroscopic and microscopic resonant tunneling to other interference 
phenomena. Consider, for instance, weak localization. If one considers a superconducting 
structure which, in the GL regime, is described by Eq. @ with a weakly disordered effective 
potential V{v), could we expect that quantum interference might cause a reduction of, i.e., 
the critical current? From the analysis presented here, the answer seems to be no, at least 
when the distance between "impurities" is much larger than ^(T): the nonlinear term would 
damp interference effects to the point of making each impurity act as isolated. In particular, 
one should not expect to observe the macroscopic equivalent of Anderson localization. A 
similar conclusion holds in principle for other mesoscopic phenomena based on geometry- 
induced quantum interference such as the modulation of current by a tuning stub:ii one 
should not expect any type of interference effect associated to the collective wave function 
in a superconductor enclosed in a specific (topologically trivial) geometry, as long as the 
relevant length scales are of order d ^ ^{T). As indicated in the Introduction, real analogs 
are to be found only for those "interference" phenomena based on fundamental symmetries 
such as gauge invariance, which gives rise to the Aharonov-Bohm effect. A qualitative 
conclusion is that care must be taken when developing physical intuitions based on the idea 
that the superconducting order parameter plays the role of a macroscopic quantum wave 
function. 
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APPENDIX A: DEDUCTION OF THE BARRIER MODEL FOR A SNS 

STRUCTURE 

Suppose we are given a quasi-one dimensional SNS structure. Let the normal metal 
ocuppy the region [— L/2,L/2]. T^^^ and TJ'^) are the critical temperatures of N and S, 
respectively. We have T]^) <T < and assume -T^^^ < T^^\ Following Refs.0@, 
we assume that both N and S have the same order parameter to gap ratio, as well 

as the same quasiparticle mass. It can be provedB that, in these conditions, the Ginzburg- 
Landau equations are valid for all x, and that the matching conditions are determined by 
the continuity of ip and its first derivative at x = ±L/2. Thus one has 
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h'^ cPR ^ o m?^ , , , / . s 

^^ + «.i? + /3.i?3 + _^.0 N>L/2 (Al) 

^ '^'^ + a^i? + /3;vi?' + ^^ = |x|<L/2 (A2) 



where = |'?/'(^)| and aiy,^^ {as,Ps) refer to the Ginzburg-Landau parameters of the 
normal (superconducting) metal, all of which depend on temperature. 

As T approaches Tj:^\ we can neglect the term jS^R^ in Eq. (A2) and the remaining 
parameters of the normal metal can be replaced by their fixed values at the T^^\ Shifting 
to reduced units of the superconductorQ we write 

-0--R + fl' + ^ = O kl>'/2 (A3) 
where I = L/^[T). This model is identical to that considered by JacobsonEl with 



ejT) 



e = > 1. (A5) 



If ^ ^ j'^/R'^, the integration of (A4) along the normal segment yields 

R'{l/2) - R'{-l/2) ^ WR, (A6) 

where \/6l = L/^n <^ 1 has been assumed, so that R{x) can be approximated as constant 
within the integral. In these conditions, the effect of the normal metal can be exactly 
mimicked by a delta barrier located at a; = with a strength 



y = to = 

as we wished to prove. 



g^W= (AT) 

SAT 



APPENDIX B: MATCHING IN THE JOSEPHSON LIMIT 

If, from the parameters intervening in Eq. (^), we define 

2 = tanh(fc(a + (Bl) 
it is not difficult to see that yo and ?/q defined in the text are given by 

yo = Z + Uz^ (B2) 



y', = V2LP{z'-z). (B3) 
So, the condition yQ < implies < z < 1. Eq. (|T^ for Seq can then be rewritten as 
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5eo{z) = I {V2lP{z' -z)+ g{Z + Uz')) . (B4) 

The additional requirement feo < implies that z is further restricted to the range Zmm < 
z < Zmax, where Zmin, Zmax the roots of be^iz) lying between and 1. We define uj as 

SO that, obviously, < < 1. 

In the Josephson limit (gf — > oo) it is easy to see that Z ~ J"^ /2g'^ [J = 2gj) and f/ ~ 1. 
Since, in the same limit, the only surviving finite root of SeqIz) is z = 0, we can neglect z^ 
in front of z in Eq. (|B4|). Then Zmin and z^nax can be easily calculated as the roots of 

- V2z + — + gz^ = (B6) 

To leading order in 1/g, we get 



•'max,min 



and, from (B5), 



1 ± Vr^J2 (B7) 



i + (2^-i)vr^ 1 

2; = j= + o(-), (B8) 

where 0(1/5'") stands for any expression such that limg_^oofi'" o{l/g"') = 0. From Eq. (B4), 
we obtain, also to leading order, 

5eoiz) = l(^V2z+^+ gz"^^ (B9) 

which can be easily shown to lead to 

Se = -6eo = uj{l - uj){l ~ J^) (BIO) 



For the other quantities in Eq. (|l^), we obtain 



J' 1 ]_ 

2{l-46e)g^^''^g^^ 

62 = 1 -2v^+o(l) 

eg = 1 + 2Vfe + 0(1) (Bll) 



The quantities /3, Qq, Qi, and m are defined as: 



Fo = ^ + o(-) (B12) 
9 9 

2a('-^Y'^no + ^ + o{h (B13) 
V 2 / 9 9 

_ 62-61 1 -2Vfe , 

m= = ——-^ + 0(1) (B14) 

63-61 1 + 2vfe 
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In this limit, the parameter a of Eq. ([Tq) becomes +1 and the global matching equation 
(pO|) can be rewritten as 

y'{-a+) + y'{a-) = g{y{-a) - y{a)). (B15) 
Solving this equation order by order in 1/(7, we get: 

sn2(l]o|m) = 0, (B16) 

which is equivalent to 



nQ = 2nK\- — n=l,2,... (B17) 
From Eqs.( |Bll| ) and ( |B13| ) we deduce 

fio ^ 2a (^^1^) (B18) 

for g large, and find that Eq. ( [B17| ) leads to Eq. (^) in the text, which we wished to prove. 

We also find that VLi can take values VLi = —2(3 (for the two symmetric solutions) or 
VLi = —Qo/a (for the two asymmetric ones). Solutions come in pairs because of the quadratic 
character of Eq. ( piO| ). The behavior under the transformation x —>■ —x can be deduced 
from the expression for the solution, which to leading order in 1/g takes the form 



yn{x) ~ ei + 62 sn^ I — ^— — — + nK + 



(2nK + VLi/g 2/5 + ^]l 



V 2a 2g 
where K stands for the elliptic integral in ( |B17| ). 



2Vfe 



(B19) 



APPENDIX C: OFFSET IN THE JOSEPHSON LIMIT 



One can insert the asymptotic expressions obtained for E?'{x) in the previous Appendix 
into the second integral of Eq. (||) and make use of the identity: 

rMK{m)-P/g du 



lini g 1 

g^QO JNK{m)+a/g 1 + g^'JSn'^ {u\m) 



,7r 



TT 



— (- - arctan(a^))5Ar + (- - aTctan{(3y/^))6M + ttN, 



(CI) 
(C2) 

7 V2 .. .. 2 .. .. y 

valid for integers A^, M and arbitrary positive values of a, (3, 7, with < m < 1. In Eq. 
Nf, is the number of even integers p satisfying N < p < M, and 6]\f is defined as one (zero) 
for N even(odd). After some lengthy but straightforward algebra, one can prove that the 
current-phase relation for the solutions of Eq. (^) is ( J > 0) 

(n + l)7r — Aip 



J{Aip) = Vl - 4fe 



sm 



Aip G [{n- l)7r, (n + l)7r] 



for the symmetric solutions, and 



A(p = niT, for all J, 



(C3) 
(C4) 



for the asymmetric solutions. 
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FIGURES 



FIG. 1. (a) Potential energy v{R) (see Eq. (13) in the text) for the mechanical analog of Eq. 
(4) with j = 0.01, g = 2 and a = 2, for A = 1 (solid line) and A = 0.8 (dotted line); the straight 
lines depict one possible trajectory of the equivalent particle that begins and ends at point B, the 
numbers in the arrows indicating time-ordered flights between kicks (the distance between the two 
horizontal lines is exactly 6e). (b) Inverted representation of the particular solution R{x) of Eq. 
(4) whose equivalent mechanical trajectory is that depicted in (a). 

FIG. 2. Critical current as function of the semidistance between barriers for several values of 
the parameter g. The distance is given in units of (,{T) and the current is given in units of the 
critical current for a single barrier with the same g (see text). 

FIG. 3. Representation of the matching equation (22) in the text. The quantity nG{5e)/TT 
is plotted for several values of the index n. A set of four solutions exists for each combination of 
n and 5e that meets the requirement nG{5e) = a, where a is the reduced semidistance between 
barriers. 

FIG. 4. Current as a function of the total phase difference for a = 1.74 (a) and a = 3.5 (b), 
in the limit of g very large. The current is given in units of the critical current for a single barrier 
with the same g, so that \im.g^^ 2gj is plotted. 

FIG. 5. Current as a function of the phase offset for several values of g in the cases a = 1.74 
(a) and a = 3.5 (b), in the same units as in Fig. 2. 

FIG. 6. Temperature dependence of the critical current in the vicinity of the critical temper- 
ature, for SNS and SNSNS structures without current concentration. d/^Q = 10 has been taken, 
and units are such that the the prefactor of (1 — T/Tc)^ in Eq. (31) equals unity. 
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